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ABSTRACT 

We  employ  complex  continuous  wavelet  transforms 
to  develop  a  consistent  mathematical  framework  capable 
of  quantifying  both  cerebrospinal  compensatory  reserve 
and  cerebrovascular  pressure-reactivity.  The  wavelet  gain 

ri'^pif)  5  defined  as  the  frequency  dependent  ratio  of  time 

averaged  wavelet  coefficients  of  intracranial  pressure 
(ICP)  and  arterial  blood  pressure  (ABP)  fluctuations, 
characterizes  the  dampening  of  spontaneous  ABP  oscilla¬ 
tions.  This  gain  is  introduced  as  a  novel  measure  of  cere¬ 
brospinal  compensatory  reserve.  For  a  group  of  12  pa¬ 
tients  who  died  as  a  result  of  cerebral  lesions  (Glasgow 

Outcome  Scale  GOS  =1)  the  average  gain  rj =  0.45 
calculated  at  y=0.05  Hz  significantly  exceeds  that  of  17 
patients  with  GOS=2:  V^abp  -0.24  with  p-  4xlCT5 
(Kruskal- Wallis  test).  We  also  study  the  dynamics  of  in¬ 
stantaneous  phase  difference  A 0(f)  between  the  fluctua¬ 
tions  of  ABP  and  ICP.  The  time-averaged  synchronization 

index  y^Bptf)  =  (sin  A^)2  +  (cos  A^)2  ,  which  depends 

upon  frequency,  yields  information  about  the  stability  of 
the  phase  difference  A</>  and  is  used  as  a  cerebrovascular 
pressure-reactivity  index.  For  both  groups  of  patients  the 
synchronization  is  strong.  We  hypothesize  that  in  patients 
who  died  the  impairment  of  cerebral  autoregulation  is 
followed  by  the  breakdown  of  residual  pressure  reactivity. 

For  example  at  f=0.05  Hz  y‘abp  =  O-70  (GOS=l)  and 

Yabp  =  0-58  (GOS=2).  While  these  two  synchronization 
levels  are  not  statistically  different  (/?=0.15)  the  corre¬ 
sponding  average  phase  difference  for  GOS=l  is  equal  to 
10  degrees  in  sharp  contrast  to  the  mean  value  of  44  de¬ 
grees  for  patients  with  GOS=2  (  p  =  1  x  1(T4  ). 


1.  INTRODUCTION 

In  the  last  decade  the  list  of  clinical  parameters  em¬ 
ployed  in  the  management  of  patients  with  severe  brain 
injuries  has  been  steadily  growing.  Continuous  measure¬ 
ment  of  ECG,  arterial  and  central  venous  pressure,  pulse 
oximetry  (Sp02),  and  end-tidal  carbon  dioxide  provide 
standard  assessment  of  systemic  cardiovascular  and  respi¬ 
ratory  functions.  These  physiologic  signals  may  now  be 
monitored  concurrently  with  intracranial  pressure  (ICP), 
jugular  venous  bulb  oxygen  saturation  (Jv02),  flow  ve¬ 
locity  in  major  cerebral  arteries,  regional  brain  tissue  con¬ 
centration  of  02,  C02,  pH,  glucose,  lactate,  and  pyruvate, 
see,  for  example,  (Kett- White  et  al.  2002)  and  references 
therein.  There  is  a  consensus  that  the  technology  used  to 
acquire  such  physiologically  diverse  time  series  has  out¬ 
paced  algorithms  used  to  analyze  them.  For  example,  con¬ 
tinuous  ICP  monitoring  has  been  an  integral  component 
of  the  management  of  patients  with  traumatic  brain  inju¬ 
ries  or  intracerebral  hemorrhages  for  over  two  decades. 
Until  recently  the  main  objective  of  such  monitoring  was 
confined  to  the  maintenance  of  adequate  cerebral  perfu¬ 
sion  pressure.  Certainly,  severe  consequences  of  even 
transient  cerebral  ischemia  to  some  extent  provide  justifi¬ 
cation  for  the  status  quo.  Nevertheless,  recent  studies 
have  clearly  demonstrated  the  clinical  significance  of 
more  sophisticated  hemodynamic  metrics  derived  from 
ICP  waveforms. 

For  example,  Czosnyka  et  al  (Czosnyka  et  al.  1996b) 
used  the  moving  correlation  coefficient  between  the  fun¬ 
damental  harmonic  of  ICP  pulse  wave  and  the  mean  ICP 
(RAP  index)  to  assess  pressure-volume  compensatory 
reserve.  In  the  same  vein,  another  ICP-derived  metric  - 
the  pressure-reactivity  index  (PRx)  (Czosnyka  et  al.  1997) 
quantifies  cerebrovascular  pressure-reactivity  with  the 
strength  of  linear  correlation  between  fluctuations  of 
arterial  blood  pressure  (ABP)  and  ICP.  A  negative  value 
of  PRx  reflects  a  normally  reactive  vascular  bed,  as  ABP 
waves  provoke  inversely  correlated  waves  in  ICP.  Posi¬ 
tive  values  of  this  measure  are  associated  with  a  passive 
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behaviour  of  a  non-reactive  vascular  bed.  This  index  was 
reported  to  correlate  with  indices  of  autoregulation  based 
on  transcranial  Doppler  ultrasonography  (Czosnyka  et  al. 
1996a;  Lang  et  al.  2002).  Furthermore,  abnormal  values 
of  both  PRx  and  RAP  indicate  poor  autoregulation  or  re¬ 
duced  cerebrospinal  compensatory  reserve,  respectively. 
These  abnormal  values  have  been  shown  to  be  predictive 
of  a  poor  outcome  following  head  injury  (Balestreri  et  al. 
2004;  Czosnyka  et  al.  1997;  Czosnyka  et  al.  1996b;  Stei¬ 
ner  et  al.  2002). 

However,  recent  studies  (Giller;  Mueller  2003;  Latka 
et  al.  2005;  Mitsis;  Marmarelis  2002;  Mitsis  et  al.  2002, 
2006;  Panerai  et  al.  1999)  have  emphasized  nonlinear, 
non-stationary  and  multiscale  aspects  of  cerebral  hemo¬ 
dynamics.  The  goal  of  this  paper  is  to  present  a  novel 
perspective  of  the  frequency-dependent,  non-stationary 
aspects  of  intracranial  pressure  dynamics.  We  adopt  and 
extend  the  mathematical  framework  of  complex  continu¬ 
ous  wavelet  transforms  to  analyze  the  dynamics  of  the 
instantaneous  phase  difference  between  the  fluctuations  of 
ABP  and  ICP  (Latka  et  al.  2005).  We  hypothesize  that  the 
frequency  dependent  strength  of  synchronization  between 
these  two  hemodynamic  time  series  reflects  pressure  reac¬ 
tivity  of  the  cerebrovascular  bed.  We  also  test  whether 
wavelet  gain  which  characterizes  the  dampening  of  spon¬ 
taneous  arterial  blood  oscillations  within  intracerebral 
space  can  quantify  cerebrospinal  compensatory  reserve. 


2.  METHODS 

The  study,  approved  by  the  Opole  Human  Research 
Committee,  comprised  37  patients  (29  males  and  8  fe¬ 
males)  treated  at  the  Department  of  Neurosurgery  of  Re¬ 
gional  Medical  Center  in  Opole,  Poland.  The  median  age 
of  subjects  was  55  ±  15  years.  In  the  cohort  19  patients 
were  diagnosed  with  subdural  haematoma,  7  with  sponta¬ 
neous  intracerebral  haemorrhage,  and  6  with  subarach¬ 
noid  hemorrhage  from  ruptured  aneurysm.  There  were  4 
cases  of  cerebral  contusions  and  a  single  case  of  traumatic 
subarachnoid  haemorrhage.  The  written  informed  con¬ 
sent  requirement  was  waived  since  only  routinely  moni¬ 
tored  variables  were  used  in  this  research. 

All  of  the  patients  were  unconscious  at  admission. 
Their  initial  neurological  condition  was  assessed  with 
Glasgow  Coma  Scale  (GCS).  5  patients  had  GCS=3;  14 
had  GCS=4;  8  had  GCS=5;  5  had  GCS=6;  4  had  GCS=7; 
and  1  had  GCS=8.  In  all  but  four  cases  the  invasive  in¬ 
tracranial  pressure  monitoring  followed  neurosurgical 
procedures.  The  tip  of  the  ICP  strain  gauge  transducer 
(Codman)  was  implanted  intraparenchymally  into  the 
frontal  lobe,  contralateraly  to  the  cerebral  lesion.  The  arte¬ 
rial  blood  pressure  was  continuously  measured  with  the 
radial  artery  cannulation  (Gabarith).  The  simultaneous 


monitoring  of  ICP,  ABP  and  ECG  was  performed  for  up 
to  7  days  after  ICP  gauge  implantation. 

Postoperative  treatment  involved  the  mechanical 
ventilation  and  optional  barbiturane  administration.  The 
cerebral  perfusion  pressure  (CPP)  was  generally  main¬ 
tained  above  70  mmHg.  The  treatment  was  scored  with 
Glasgow  Outcome  Scale  (GOS).  12  out  of  37  patients 
died  (GOS=l).  For  the  other  patients  the  scores  were  dis¬ 
tributed  as  follows:  17  with  GOS=2;  6  with  GOS=3;  1 
with  GOS=4;  and  1  with  GOS=5.  The  statistical  analysis 
was  confined  to  the  two  largest  subgroups,  i.e.,  those  with 
GOS=l  and  those  with  GOS=2.  The  physiologic  data 
acquired  during  monitoring  of  patients  with  better  clinical 
condition  (GOS  >  2)  were  used  for  case  studies. 


Beat-to-beat  average  values  of  blood  and  intracranial 
pressures  were  calculated  via  waveform  integration  of  the 
corresponding  signals  sampled  at  100  Hz  and  digitized  at 
16  bits.  In  numerical  calculations  non-uniformly  spaced 
time  series  were  resampled  at  2  Hz  with  the  help  of  cubic 
spline  interpolation. 


The  wavelet  transform  is  an  integral  transform  which 
employs  basis  functions,  known  as  wavelets,  localized  in 
both  time  and  frequency.  Such  wavelets  are  constructed 
from  a  single  mother  wavelet  y/(t)  by  means  of  transla¬ 
tions  and  dilations.  The  wavelet  coefficient  Ws(a,t0)  is 
obtained  by  filtering  the  signal: 

OO  ✓  \ 

ws(a,t0)=-Ljv,*(^J.S(t)dt  (i) 

— OO 

and  quantifies  the  similarity  of  the  signal  s(t )  to  the 


wavelet  dilated  by  scale  a  and  centered  at  t0  (the  asterix 
denotes  complex  conjugation).  The  instantaneous  phase 
(/){t)  associated  with  the  signal  s(t)  can  be  readily  ex¬ 


tracted  by  calculating  its  wavelet  transform  with  a  com¬ 
plex  mother  function  to  obtain 


</>(a,t0)  =  -/log 


KM 

\wsM\ 


(2) 


The  dual  localization  of  wavelets  makes  it  possible  to 
associate  a  pseudofrequency  fa  with  the  scale  a  : 


/«  = 


fc_ 

aSt 


(3) 


where  fc  is  the  center  frequency  and  St  is  the  sampling 
period.  Herein,  we  employ  the  complex  Morlet  wavelet: 


v(t)  = 


1  c2  nifcte-?!fb 

^/7 Tb 


(4) 


and  set  both  the  center  frequency  fc  and  the  bandwith 


parameter  fb  to  1.  The  detailed  description  of  the  nu¬ 
merical  algorithms  used  in  wavelet  analyses  may  be 
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found  in  (Latka  et  al.  2003;  Latka  et  al.  2005)  and  will  not 
be  reproduced  here. 


We  quantify  the  interdependence  of  ABP  and  ICP 
time  series  fluctuations  with  the  scale  (frequency)  de¬ 
pendent  synchronization  index: 

?ABp(a) =  (sin  A(2>)2  +(cos  A^)2  (5) 

where  A</>  =  (/)ICP  -  (j)ABP  is  the  instantaneous  phase  dif¬ 
ference  between  the  ABP  and  ICP  time  series.  The  angu¬ 
lar  brackets  denote  averaging  over  time.  The  synchroniza¬ 
tion  index  lies  in  the  interval  0  <  y  <  1 .  A  vanishing  index 

y  -  0  corresponds  to  a  uniform  distribution  of  the  phase 


differences  (no  synchronization)  while  y  -  1  corresponds 

to  perfect  synchronization  (phase  locking  of  the  two  proc¬ 
esses). 

In  the  most  straightforward  approach  we  define  the 
scale-dependent  time-averaged  wavelet  gain  rj  as 

(K/cpOMo)  I) 


Vabp (fl) — 


l\WABP{a,t0)\) 


(6) 


time  [s] 


Fig.  1.  Wavelet  analysis  of  arterial  (left  column)  and  intrac¬ 
ranial  (right  column)  pressure  fluctuations  for  a  patient  with 
mild  hypertension.  First  row:  The  spontaneous  fluctuations 
averaged  over  a  cardiac  beat.  Second  row:  time  evolution  of  the 
modulus  of  the  wavelet  coefficient  for  scale  a=10  (this  scale 
corresponds  to  pseudofrequency  0.2  Hz).  Third  row:  the  den¬ 
sity  plot  of  the  modulus  of  the  wavelet  coefficients.  The  hori¬ 
zontal  color  bar  at  the  bottom  of  the  figure  reveals  the  assign¬ 
ments  of  colors  to  values  of  wavelets  coefficients.  All  coeffi¬ 
cients  with  modulus  greater  than  4  are  marked  as  dark  brown. 
Strong  damping  of  ABP  fluctuations  is  characteristic  of  physio¬ 
logical  cerebral  hemodynamics. 


The  gain  quantifies  the  damping  of  spontaneous  ABP 
fluctuations.  We  hypothesize  that  the  synchronization 
index  y  is  determined  by  reactivity  of  the  vascular  bed 
while  the  gain  77  may  be  used  as  a  novel  measure  of  cere¬ 
brospinal  compensatory  reserve. 


3.  DISCUSSION 

It  is  apparent  that  the  fundamental  Fourier  component 
of  intracranial  pressure  waveform  has  a  frequency  equal 
to  the  heart  rate.  Czosnyka  et  al  (Czosnyka  et  al.  1996b) 
argue  that  linear  correlation  between  the  amplitude  of  the 
pulse  harmonic  (AMP)  and  mean  intracranial  pressure 
provides  information  about  accessible  intracerebral 
space.  More  specifically,  the  relation  between  ICP  and 
this  volume  is  given  by  a  nonlinear,  sigmoid  pressure- 
volume  curve.  At  low  ICP  this  curve  is  flat  indicating 
good  compensatory  reserve.  With  increasing  pressure  the 
relation  becomes  exponential.  The  second  flat  zone,  seen 


Fig.  2.  Wavelet  analysis  of  arterial  (left  column)  and  intracra¬ 
nial  (right  column)  pressure  fluctuations  for  a  patient  with  se¬ 
vere  hypertension.  First  row:  The  spontaneous  fluctuations 
averaged  over  a  cardiac  beat.  Second  row:  time  evolution  of  the 
modulus  of  the  wavelet  coefficient  for  scale  a=10  (this  scale 
corresponds  to  pseudofrequency  0.2  Hz).  Third  row:  the  den¬ 
sity  plot  of  the  modulus  of  the  wavelet  coefficients.  The  hori¬ 
zontal  color  bar  at  the  bottom  of  the  figure  reveals  the  assign¬ 
ments  of  colors  to  values  of  wavelets  coefficients.  All  coeffi¬ 
cients  with  modulus  greater  than  4  are  marked  as  dark  brown. 
Weak  damping  of  ABP  fluctuations  stems  from  severe  cerebral 
swelling 


3 


Fig.  3.  Frequency  dependence  of  wavelet  gain  for:  a  patient 
with  mild  (dashed  line)  and  severe  (solid  line)  cerebral  hyper¬ 
tension. 


at  very  high  ICP  reflects  derangement  of  normal  cere¬ 
brovascular  responses.  More  specifically,  they  define  the 
RAP  index  as  the  coefficient  of  correlation  between  40 
consecutive,  time-averaged  data  points  of  AMP  and  ICP 
acquired  over  6.4  s.  The  values  of  RAP  close  to  0  indicate 
good  pressure-volume  compensatory  reserve  at  low  ICP. 
Even  though  the  simplicity  and  clinical  relevance  of  the 
discussed  approach  is  appealing,  the  fundamental  ques¬ 
tion  arises  as  to  whether  the  response  to  changes  in  intrac¬ 
ranial  blood  volume  is  frequency  dependent. 

In  Fig.  1  and  Fig.  2  we  compare  time  evolution  of  the 
modulus  of  wavelet  coefficients  of  ABP  and  ICP  fluctua¬ 
tions  for  two  patients  with  mild  and  severe  cerebral  hy¬ 
pertension  respectively.  On  the  basis  of  these  case  studies 
one  can  hypothesize  that  strong  damping  of  spontaneous 
ABP  fluctuations  is  characteristic  of  physiologic  hemody¬ 
namics.  Thus,  the  wavelet  gain  77  which  quantifies  such 
attenuation  is  a  plausible  metric  for  pressure-volume 
compensatory  reserve. 


f  [Hz] 

Fig.  4.  Group  averaged  wavelet  gain  as  a  function  of  frequency. 
The  thick  solid  line  corresponds  to  the  patients  who  died 
(GOS=l)  and  the  dashed  line  to  those  with  GOS=2.  The  thin 
lines  represent  the  corresponding  standard  error  of  the  mean. 

In  Fig.  3  we  plot  the  gain  as  a  function  of  frequency 
for  the  cases  analyzed  in  Figs.  1  and  2.  The  elevated  value 
of  wavelet  gain  over  the  entire  frequency  range  is  evident 
for  a  subject  with  severe  cerebral  swelling.  It  is  also  worth 
pointing  out  that  the  damping  of  ABP  fluctuations  is  sig¬ 
nificantly  weaker  (as  indicated  by  higher  values  of  gain) 
in  low-frequency  part  of  the  spectrum.  These  conclusions 
are  corroborated  by  the  statistical  analysis  of  the  wavelet 
gain  for  two  group  of  patients:  those  with  GOS=l  and 
those  with  GOS=2  (cf.  Fig.  4).  With  decreasing  fre¬ 
quency,  the  mean  gain  for  both  groups  starts  to  diverge  at 
approximately  0.2  Hz.  For  example,  for  patients  who  died 

the  average  gain  77 ^  =  0.45  calculated  aty=0.05  Hz  sig¬ 
nificantly  exceeds  that  of  patients  with  GOS=2: 
Vabp  ~  0-24  ( p  =  4  x  1(T5  obtained  from  the  Kruskal- 

Wallis  test). 
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Fig.  5.  Time  evolution  of  frequency  dependent  normalized  phase  difference  (A/2  n  )  between  fluctuations  of  arterial  and  intracranial  pressure 
fluctuations  for  a  patient  with:  (a)  mild  hypertension  (cf.  Fig.  1)  and  (b)  severe  hypertension  (cf.  Fig.  2). 
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f  [Hi] 

Fig.  6.  Frequency  dependence  of  synchronization  parameter  for: 
a  patient  with  mild  (dashed  line)  and  severe  (solid  line)  hyper¬ 
tension. 


In  the  phenomenological  description  of  cerebral  he¬ 
modynamics,  fluctuations  of  ABP,  due  to  pressure  reac¬ 
tivity  of  cerebral  vessels,  lead  to  fluctuations  of  ICP.  The 
vascular  origin  of  the  ICP  slow  waves  through  intracranial 
blood  volume  variations  is  commonly  accepted 
(Steinmeier  et  al.  1996).  Nevertheless,  the  quantitative 
description  of  this  clinically  significant  effect  is  the  sub¬ 
ject  of  ongoing  research.  To  shed  new  light  on  the  inter¬ 
play  of  ABP  and  ICP  we  chose  to  employ  complex  wave¬ 
let  transforms  to  study  the  time  evolution  of  the  instanta¬ 
neous  phase  difference  between  the  two  time  series.  Fig.  5 
exemplifies  such  an  approach  for  two  previously  dis¬ 
cussed  case  studies.  Upon  inspection  of  both  density  plots 
in  this  figure  one  finds  that  the  presence  of  the  horizontal 
multicolor  strip  for  a  subject  with  mild  cerebral  hyperten¬ 
sion  (Fig.  5a)  is  just  a  manifestation  of  low-frequency 
phase  variability.  Such  variability  is  almost  entirely  lost 
during  severe  hypertension  (Fig.  5b).  Synchronization 
strength  between  fluctuations  of  ABP  and  ICP  can  be 
quantified  with  the  help  of  synchronization  parameter  y 
which  in  Fig.  6  is  plotted  as  a  function  of  frequency  for 
both  subjects.  The  very  high  values  of  y  for  the  patient 
with  GOS=l  reflect  pathological  entrainment  of  ABP  and 
ICP  time  series.  Such  entrainment  can  also  be  observed  in 
Fig.  7  which  shows  the  synchronization  parameter  aver¬ 
aged  over  both  groups.  For  this  hemodynamic  metric  the 
differences  between  the  groups  are  not  as  strongly  pro¬ 
nounced  across  the  displayed  frequency  range  as  were  the 
differences  in  wavelet  gain  (Fig.  4).  Aty^0.05  Hz  we  find 

rj™p  =  0.70  (GOS=l)  and  =  0.58  (GOS=2).  These 

two  synchronization  levels  are  not  statistically  different 
(77=0.15).  In  principle  one  might  argue  that  for  both 
groups  cerebrovascular  reactivity  was  compromised. 
However,  strong  phase-locking  allows  us  to  characterize 
cerebral  hemodynamics  also  in  terms  of  the  time- 
averaged  value  of  the  phase  difference  between  the  ABP 
and  ICP  waveforms. 


1  [Hi] 

Fig.  7.  Group  averaged  synchronization  parameter  as  a  function 
of  frequency.  The  thick  solid  line  corresponds  to  the  patients 
who  died  (GOS=l)  and  the  dashed  line  to  those  with  GOS=2. 
The  thin  lines  represent  the  corresponding  standard  error  of  the 
mean. 


Fig.  8  shows  the  group  averaged  phase  difference  as 
a  function  of  frequency.  At  f=0.05  Hz  the  corresponding 
difference  for  GOS=l  is  equal  to  10  degrees  in  sharp 
contrast  to  the  mean  value  of  44  degrees  for  patients  with 
GOS=2  (circular  variant  of  the  t-test  (Batschelet  1981) 

yielded  the  value  of  p  =  lxl(T4).  The  low-frequency 

variability  of  the  phase  difference  between  the  spontane¬ 
ous  oscillations  of  ABP  and  blood  flow  velocity  in  major 
cerebral  arteries  is  the  hallmark  of  intact  cerebral  auto¬ 
regulation  (Giller  1990;  Latka  et  al.  2005).  Therefore,  we 
hypothesize  that  in  patients  who  died  the  impairment  of 
cerebral  autoregulation  is  followed  by  the  breakdown  of 
residual  pressure  reactivity.  It  is  worth  emphasizing  that 
such  insight  is  unique  to  the  wavelet  analysis.  In  the  PRx 
approach  the  gradual  loss  of  cerebrovascular  bed  pressure 
reactivity  is  manifested  solely  by  the  index’s  value  ap¬ 
proaching  one. 


Fig.  8.  The  time  averaged  angle  between  ABP  and  ICP  fluctua¬ 
tions  for:  patients  with  GOS=l  (thick  solid  line)  and  GOS=2 
(thick  dashed  line).  The  thin  lines  delineate  95%  confidence 
intervals. 

In  the  pioneering  studies  of  cerebrospinal  fluid 
system  Marmarou  (Marmarou  1973)  postulated  an  expo¬ 
nential  volume-pressure  relation.  Such  interdependence  of 
variables  underlies  the  PVI  index  which  quantifies  the 
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compliance  over  the  whole  physiological  range  of  ICP. 
Miller  et  al  (Miller  et  al.  1973)  advocated  the  volume 
pressure  response  (VPR)  as  a  measure  of  elastance  (in¬ 
verse  of  compliance)  which  is  independent  of  functional 
relation  between  craniospinal  volume  and  pressure.  In 
both  approaches,  now  commonly  accepted  as  the  gold 
standards,  the  corresponding  indices  are  calculated  from 
the  pressure  change  resulting  from  volume  injection  or 
withdrawal  of  fluid  from  the  CSF  space.  An  exhaustive 
list  of  relevant  references  is  given  in  (Reilly,  Bullock 
2005). 

Physiological  causes  of  elevated  ICP  may  be  broadly 
classified  into  vascular  and  non- vascular.  Vascular  me¬ 
chanisms  are  associated  with  either  passive  or  activite 
vasolidation  of  the  cerebrovascular  bed.  Active  vasolida- 
tion  may  be  brought  about  by  increased  concentration  of 
C02  and/or,  within  the  operating  range  of  pressure  auto¬ 
regulation,  decreased  cerebral  perfusion  pressure.  The  list 
of  non-vascular  causes  is  long  and  includes  increased: 
brain  water  content  (edema),  intracerebral,  extradural  or 
subdural  mass  as  well  as  CSF  outflow  resistance  (Reilly; 
Bullock  2005). 

Marmarou  et  al  demonstated  that  in  patients  with 
severe  head-injury  elevated  ICP  can  to  a  large  extent  be 
attributed  to  vascular  mechanisms  (Marmarou  et  al. 
1987).  Gray  and  Rosner  determined  that  PVI  and  CPP 
may  be  linearly  anticorrelated  or  correlated  depending  on 
whether  the  CPP  range  is  beyond  or  within  the  autoregu¬ 
lating  range  (Gray;  Rosner  1987a,  1987b).  Detrimental 
influence  of  low  CPP  on  intracranial  compliance  was  re¬ 
cently  corroborated  by  Portella  et  al  (Portella  et  al.  2005). 
Interestingly  enough,  cerebrovascular  compensatory  re¬ 
serve  is  also  time-dependent.  Anile,  Portnoy,  and  Branch 
found  that  VPR  indices  derived  from  slow,  medium  and 
rapid  bolus  injection  differ  (Anile  et  al.  1987).  The  dif- 
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ferences  result  from  the  existence  of  two  distinct  com¬ 
ponents  of  craniospinal  compliance.  The  short-time  phys¬ 
ical  compliance  is  determined  by  such  factors  as  spinal 
dura  matter  expansion  while  the  time  constant  of  physio¬ 
logical  compliance  associated  with  cerebrovascular  altera¬ 
tions  (e.g.  venous  outflow  resistance)  is  significantly 
longer.  Thus,  to  fully  understand  the  intricate  interplay  of 
craniospinal  volume  and  pressure,  the  dynamic  and  the 
viscoelastic  properties  of  CSF,  nervous  tissue  and  vascu¬ 
lar  factors  must  all  be  considered. 

To  address  this  issue  in  the  most  straightforward  way, 
in  Fig.  9  the  values  of  wavelet  gain  and  synchronization 
index,  averaged  over  the  low  frequency  part  of  the  spec¬ 
trum,  were  used  as  the  coordinates  of  points  representing 
the  overall  hemodynamic  condition.  It  comes  as  no  sur¬ 
prise  that  the  high-synchronization  high-gain  region  of  the 
cerebral  hemodynamic  state  space  to  a  large  extent  coin¬ 
cides  with  high  values  of  intracranial  pressure.  Neverthe¬ 
less,  one  can  easily  find  cases  with  relatively  high  wavelet 
gain,  moderate  ICP  and  largely  intact  pressure  reactivity 
of  the  vascular  bed.  Thus,  further  studies  should  explore 
whether  the  onset  of  severe  cerebral  hypertension  is  pre¬ 
ceded  by  high-gain  episodes.  If  such  episodes  are  suffi¬ 
ciently  long  they  could  become  the  cornerstone  of  a  ICP 
hypertension  prediction  algorithm.  It  also  worth  pointing 
out  that  overnight  ICP  monitoring  may  be  used  in  the 
diagnosis  of  patients  with  hydrocephalus  symptoms  or 
shunt  malfunctions  (Czosnyka  et  al.  2004).  The  frequency 
dependent  aspects  of  cerebrospinal  fluid  dynamics  dis¬ 
cussed  herein  give  strong  indications  that  such  monitor¬ 
ing  can  provide  more  detailed  assessment  of  patients’ 
clinical  conditions  than  do  limited  time  infusion  tests. 
From  the  clinical  point  of  view,  the  studies  of  time  evolu¬ 
tion  of  the  wavelet  gain  in  patients  with  normotensive 
hydrocephalus,  are  particularly  important. 


GOS  2 
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Fig.  9.  The  values  of  wavelet  gain  and  synchronization  index,  averaged  over  the  low  frequency  part  of  the  spectrum  (0.02-0.07  Hz),  were 
used  as  the  coordinates  of  points  representing  the  hemodynamic  condition  of  patients  with  GOS=l  (left  plot)  and  GOS=2  (right  plot). 
The  vertical  color  bar  shows  the  assignment  of  colors  to  values  of  intracranial  pressure.  Each  dot  corresponds  to  15  minute  monitoring 
interval. 
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4.  CONCLUSIONS 

Diverse  approaches  have  been  used  to  elucidate  inter¬ 
relations  between  intracranial  pressure  and  other  hemo¬ 
dynamic  variables.  Most  of  the  research  effort  has  been 
focused  on  propagation  of  the  pulse  wave  from  the  arte¬ 
rial  bed  to  the  intracranial  space  (Aboy  et  al.  2003;  Ave- 
zaat  et  al.  1979;  Daley  et  al.  1995;  Pettorossi  et  al.  1978; 
Portnoy;  Chopp  1981;  Portnoy  et  al.  1982).  However, 
during  the  last  decade  the  clinical  significance  of  low  fre¬ 
quency  ICP  waves  has  been  repeatedly  emphasized 
(Steinmeier  et  al.  1996).  Herein  we  presented  the  mathe¬ 
matical  framework  capable  of  quantifying  intrinsically 
non-stationary  properties  of  cerebral  hemodynamics 
across  the  broad  range  of  physiologically  relevant  fre¬ 
quencies.  We  are  convinced  that  optimal  therapy  man¬ 
agement  of  patients  with  cerebral  lesions  hinges  on  thor¬ 
ough,  qualitative  assessment  of  their  clinical  condition. 
The  further  validation  of  the  proposed  assessment  of 
compensatory  reserve  and  pressure  reactivity  will  involve 
the  comparison  with  the  gold  standards  (Diehl  2002; 
Reilly;  Bullock  2005)  as  well  as  the  methods  which  fairly 
recently  gained  clinical  acceptance  (Piper  et  al.  1999; 
Piper  et  al.  1990).  This  issue  is  the  subject  of  our  ongoing 
research. 
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